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interacting two-dimensional electron systems are reviewed and critically analyzed. Special attention 
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Two-dimensional (2D) electron systems are realized 
when the electrons are free to move in a plane but 
their motion perpendicular to the plane is quantized in 
a confining potential well. At low electron densities in 
such systems, the strongly-interacting limit is reached be- 
cause the kinetic energy is overwhelmed by the energy of 
electron-electron interactions. The interaction strength 
is characterized by the ratio between the Coulomb en- 
ergy and the Fermi energy, r* = E ee /Ep. Assuming 
that the effective electron mass is equal to the band mass, 
the interaction parameter r* in the single-valley case re- 
duces to the Wigner-Seitz radius, r s — l/(7rn s ) 1 / 2 as and 
therefore increases as the electron density, n s , decreases 
(here ag is the Bohr radius in semiconductor). Possi- 
ble candidates for the ground state of the system include 
Wigner crystal characterized by spatial and spin order- 
ing [lj , ferromagnetic Fermi liquid with spontaneous spin 
ordering Q, paramagnetic Fermi liquid Q, etc. In the 
strongly-interacting limit (r s 3> 1), no analytical theory 
has been developed to date. According to numeric sim- 
ulations Q, Wigner crystallization is expected in a very 
dilute regime, when r s reaches approximately 35. The 
refined numeric simulations Q have predicted that prior 
to the crystallization, in the range of the interaction pa- 
rameter 25 < r s < 35, the ground state of the system 
is a strongly correlated ferromagnetic Fermi liquid. At 
higher electron densities, r s ~ 1, the electron liquid is 
expected to be paramagnetic, with the effective mass, m, 
and Lande g factor renormalized by interactions. Apart 
from the ferromagnetic Fermi liquid, other intermediate 
phases between the Wigner crystal and the paramagnetic 
Fermi liquid may also exist. 

In real 2D electron systems, the inherent disorder leads 
to a drastic change of the above picture, which signifi- 
cantly complicates the problem. According to the scaling 
theory of localization |y] , all electrons in a disordered in- 
finite noninteracting 2D system become localized at zero 
temperature and zero magnetic field. At finite tempera- 
tures, regimes of strong and weak localizations are distin- 
guished: (i) if the conductivity of the 2D electron layer is 
activated, the resistivity diverges exponentially as T — > 0; 
and (ii) in the opposite limit of weak localization the 
resistivity increases logarithmically with decreasing tem- 
perature — an effect originating from the increased prob- 
ability of electron scattering from impurities back to the 



2 



0) 
"O 

O , 
03 



WC 




.'FFL' 

, Paramagnetic 

' Fermi liquid 



i / 

_1 !_ 



Electron density 



FIG. 1: Schematic phase diagram in disorder vs n s plane. 
The Wigner crystal (WC) regime is preceded by the ferro- 
magnetic Fermi liquid (FFL) p]. The high-disorder region is 
shaded. 



starting point. Interestingly, the incorporation of weak 
interactions (r s < 1) between the electrons promotes the 
localization 0. However, for weak disorder and r e >l a 
possible metallic ground state was predicted |SL ISL Il0| - 

In view of the competition between the interactions 
and disorder, high- and low-disorder limits can be consid- 
ered. Since many of the experimental groups have made 
little distinction between these, there has been confusion 
about interpretation of the data. In highly-disordered 
electron systems, the range of low densities is not ac- 
cessible as the strong (Anderson) localization sets in (see 
Fig.nj. At higher densities, a logarithmic- in-T correction 
to the resistivity was observed in numerous experiments 
(see, e.g., Refs. 0, 0, El), providing support for the 
weak localization theory. Apparently, extrapolation of 
the weak corrections to T = is not justified and, there- 
fore, those studies cannot serve as confirmation of the 
scaling theory. The latter is remarkable by the princi- 
pal impossibility of experimental verification because all 
experiments are performed in samples with finite dimen- 
sions at finite temperatures. The question whether or not 
the scaling theory works is essentially a matter of belief. 

The case of low-disordered electron systems is much 
more interesting. Low electron densities corresponding 
to the strongly-interacting limit become accessible. Ex- 
perimental results on the metal-insulator phase diagram 
in perpendicular magnetic fields revealed a close simi- 
larity between the insulating phase at low densities and 
the quantum Hall states. Thus, they excluded the for- 
mation of a pinned Wigner crystal in available samples, 
but supported the existence of a metallic state in zero 
field [3 E E3. As the magnetic field is decreased, 
the extended states in the Landau levels were observed 
to float up in energy relative to the Landau level cen- 
ters and merge to form a metallic state in the B = 
limit. This contradicts the theoretical scenario that in 
the limit of zero magnetic field the extended states should 
float up indefinitely in energy [TtI ITsj leading to an in- 
sulating ground state. The metallic state was found to 



be remarkable by the strong drop of resistivity with de- 
creasing temperature 0, |2£l . Although the origin of 
the effect was attributed to strong electron-electron in- 
teractions, the underlying physics remained unclear until 
recently. 

A breakthrough in understanding this topic occurred 
in the past four years. After a strongly enhanced ratio 
gm of the spin and the cyclotron splittings was found 
at low n s in silicon MOSFETs [22j. it became clear 
that the system behaves well beyond the weakly inter- 
acting Fermi liquid. It was reported that the mag- 
netic field required to produce complete spin polariza- 
tion, B c oc n s /gm, tends to vanish at a finite electron 
density n x w 8 x 10 10 cm~ 2 , which is close to the criti- 
cal density, n c , for the metal-insulator transition in this 
electron system [23|, H3, H3- These findings point to a 
sharp increase of the spin susceptibility, \ 00 S m i an d 
possible ferromagnetic instability in dilute silicon MOS- 
FETs. In very dilute GaAs/AlGaAs heterostructures, a 
similar behavior has been observed in both 2D hole and 
2D electron systems [2(| [SIJ . Recently, experimental re- 
sults have indicated that in silicon MOSFETs it is the 
effective mass, rather than the g factor, that sharply in- 
creases at low electron densities (28|. They have also 
indicated that the anomalous rise of the resistivity with 
temperature is related to the increased mass. The mag- 
nitude of the mass does not depend on the degree of spin 
polarization, which points to a spin-independent origin of 
the effective mass enhancement [2^, H(J • It is interesting 
that the observed effects are more pronounced in silicon 
MOSFETs compared to GaAs/AlGaAs heterostructures, 
although the fractional quantum Hall effect, which is usu- 
ally attributed to electron-electron interactions, has not 
been reliably established in silicon MOSFETs. 

The fact that n x is close to the critical density n c indi- 
cates that the metal-insulator transition in silicon sam- 
ples with very low disorder potential is a pro pert y of a 
clean 2D system and is driven by interactions [23] . This 
is qualitatively different from a localization-driven tran- 
sition in more-disordered samples that occurs at appre- 
ciably higher densities than n x which are also dependent 
on disorder strength. In this review, attention is focused 
on results obtained in the clean regime. 

There are several good reviews on the topic in ques- 
tion (see, e.g., Refs. 0, H2). However, they have ei- 
ther become outdated or do not express criticism toward 
the items of general belief, e.g., the scaling theory. Be- 
low, I describe the main experimental results and give a 
broad view on the metal-insulator transition and anoma- 
lous properties of 2D electron systems at low densities. 



II. METAL-INSULATOR PHASE DIAGRAMS 
IN A MAGNETIC FIELD 

Metal-insulator transitions in perpendicular magnetic 
fields attracted a great deal of interest in the past 
decade. The experimental activity was strongly stim- 
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ulated by theoretical predictions that Wigner crystal- 
lization is promoted in the presence of a magnetic field 
[IE EH EE HE llll ■ The insulating phase at low elec- 
tron densities was mainly studied whose origin was at- 
tributed to possible formation of the Wigner crystal 
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However, the latter has been precluded in studies of the 
metal-insulator phase diagram including quantum Hall 
states, which show a close similarity of all insulating 
phases in available samples fl4L Il5l flr| . It is interest- 
ing that there are some firmly-established experimental 
results which have not attracted much of the theorists' 
attention. These include (i) oscillations of the metal- 
insulator phase boundary as a function of perpendicular 
magnetic field; and (ii) finite bandwidth of the extended 
states in the Landau levels. 




A. Floating-up of the extended states in 
perpendicular magnetic fields 

1. First observation 

The scaling theory of localization was challenged by the 
quantum Hall effect (quantization of the Hall resistivity, 
Pxy = h/ve 2 , at integer filling factor v accompanied with 
vanishing longitudinal resistivity, p xx ) [5j| which implies 
the existence of extended states in the Landau levels (see 
section Hi C|) . To reconcile these two, it was theoretically 
predicted almost immediately that the extended states in 
the Landau levels cannot disappear discontinuously with 
decreasing magnetic field but must float up indefinitely in 
energy in the limit of B = 0, 0] . The expected phase 
diagram is shown in the inset to Fig.[2Ia). An equivalent 
diagram plotted in disorder versus inverse filling factor 
{l/v = eB/hcn s ) plane is known as the global phase 
diagram for the quantum Hall effect j54[. As long as no 
merging of the extended states was considered to occur, 
their piercing of the Fermi level was predicted to cause 
quantization of the Hall conductivity in weak magnetic 
fields miHI. 

The first attempt [fif to experimentally determine the 
metal-insulator phase diagram at low temperatures in 
low-disordered silicon MOSFETs already revealed dis- 
crepancies with the theory (see Fig.^a)). In that paper, 
a somewhat arbitrary criterion for the longitudinal con- 
ductivity, a xx — e 2 /20ft., was used to map out the phase 
boundary that corresponds to the Anderson transition 
to the regime of strong localization. However, firstly, 
the phase boundary was shown to be insensitive to the 
choice of the cutoff value (see, e.g., Ref. [13]), and sec- 
ondly, that particular cutoff value is consistent with the 
results obtained for quantum Hall states by a vanishing 
activation energy combined with a vanishing nonlinearity 
of current-voltage characteristics when extrapolated from 
the insulating phase ^| (note that for the lowest-density 
phase boundary, a lower value a xx ~ 100 kOhm at a tem- 
perature « 25 mK follows from the latter method). The 
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FIG. 2: (a) Metal-insulator phase diagram in a low- 
disordered 2D electron system in silicon MOSFETs obtained 



using a cutoff criterion, a x 



e /20h, at a temperature 



~ 25 mK. The dimensionless a xy h/e 2 in different insulating 
phases is indicated. The slope of the dashed line is close 
to e/2hc. A sketch of the expected phase diagram is dis- 
played in the inset. From Ref. 14]. (b) The map of extended 
states determined by maxima in a xx in a low-disordered sili- 
con MOSFET. Numbers show a xy in units of e 2 /h. Adopted 
from Ref. El. 



metallic phase surrounds each insulating phase as charac- 
terized by the dimensionless Hall conductivity, a xy h/e 2 , 
that counts the number of quantum levels below the 
Fermi level (in bivalley (1 00)-silicon MOSFETs, spin and 
valley degeneracies of the Landau level should be taken 
into account). This indicates that the extended states 
indeed do not disappear discontinuously. However, with 
decreasing magnetic field they float up in energy relative 
to the Landau level centers and merge forming a metallic 
state in the limit of B = (see sections IIIBI and III D|l . 
Besides, the phase boundary at low electron densities os- 
cillates as a function of B with minima corresponding 
to integer filling factors. The phase boundary oscilla- 
tions manifest themselves in that at electron densities 
near the B = metal-insulator transition, the magne- 
toresistance oscillates with an amplitude that diverges as 
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T — > 38] ; the regions in which the magnetoresistance 
diverges are referred to as reentrant insulating phase (see 
section Hi B|) . 



2. Other methods and 2D carrier systems 

The topology of the observed metal-insulator phase 
diagram — merging the extended states and, hence, 
the presence of direct transitions between the insulat- 
ing phase with a xy = and quantum Hall phases with 
a xy h/e 2 > 1 — is robust being insensitive to the method 
for spotting the phase boundary and to the choice 

of 2D carrier system |Si||6(|. It was verified using a crite- 
rion of vanishing activation energy and vanishing nonlin- 
earity of current-voltage characteristics as extrapolated 
from the insulating phase, which allows more accurate de- 
termination of the Anderson transition • Kravchenko 
et al. also applied for similar silicon MOSFETs a 
method that had been suggested in Ref. |6lJ . They stud- 
ied extended states by tracing maxima in the longitudi- 
nal conductivity in the (B,n s ) plane (see Fig.^b)) and 
found a good agreement with the aforementioned results. 
A similar merging of at least the two lowest extended 
states was observed in more-disordered 2D hole systems 
in a GaAs/AlGaAs heterostructure (see Fig. Eta)) 
and in a Ge/SiGe quantum well |(5(j (see Fig. EJb)). In 
the former case the extended states were determined by 
peaks in a xx or temperature-independent crossing points 
in p xx ; in the latter they were associated either with 
maxima in p xx and/or dp xy /dB, or with crossing points 
of p xx at different temperatures. It is noteworthy that 
a bad combination of the criterion for determining the 
phase boundary and the 2D carrier system under study 
may lead to a failure of mapping out the phase diagram 
down to relatively weak magnetic fields. In Ref. [6l| . 
extended states were studied by measuring maxima in 
the longitudinal conductivity in the (£?, n s ) plane for the 
strongly-disordered 2D electron system in GaAs/AlGaAs 
heterostructures (see Fig.^c)). Because of strong damp- 
ing of the Shubnikov-de Haas oscillations in low magnetic 
fields, the desired region on the phase diagram below 2 T 
was not accessible in that experiment. This invalidates 
the claim of Glozman et al. [61| that the extended states 
do not merge . The behavior of the lowest extended 
state in Fig.[3Jc), which Glozman et al. [6l| claim to float 
up above the Fermi level as B — ► 0, simply reflects the 
occurrence of a phase boundary oscillation minimum at 
filling factor v — 2, similar both to the minimum at v = 1 
in Fig. Ola) and to the case of silicon MOSFETs (Fig.EJ). 
Such a minimum manifests itself in that there exists a 
minimum in p xx at integer v > 1 that is straddled by 
the insulating phase [H Hill El S E3 • To this end, 
all available data for the metal-insulator phase diagrams 
agree well with each other, except those in the vicinity 
of B = 0. In weak magnetic fields, experimental results 
obtained in 2D electron systems with high disorder are 
not method- independent. Glozman et al. |61| found that 
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FIG. 3: (a) A map of the extended states for a highly- 
disordered 2D hole system in GaAs/AlGaAs heterostructures. 
Each data point represents a distinct peak in a xx or a tem- 
perature independent crossing point in p xx . Numbers show 



the value of a xy h/e . Adopted from Ref. 



(b) A map of 



the extended states for a highly-disordered 2D hole system in 
a Ge/SiGe quantum well. The open circles represent maxima 
in p xx and/or dp xy /dB. The solid circles correspond to cross- 
ing points of p xx at different temperatures. Numbers show 
the value of a xy h/e 2 . Adopted from Ref. |g(J. (c) Behav- 
ior of the extended states determined by maxima in a xx in a 
strongly-disordered 2D electron system in GaAs/AlGaAs het- 
erostructures. Numbers show a xy in units of e 2 jh. Adopted 
from Ref. EJ. 



the cutoff criterion yields basically a flat phase boundary 
towards B = 0, which is in agreement with the data for 
silicon MOSFETs (Fig. Ufa)). On the contrary, Hilke et 
al. |6d ] employed the method based on temperature de- 
pendences of p xx and obtained a turn up on the phase 
boundary in Fig. Otb) . Note that the validity of the data 
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for the lowest extended state at magnetic fields < 1.5 T in 
Fig. EI b) is questionable because the weak temperature 
dependences of p xx as analyzed by Hilke et al. 60] cannot 
be related to cither an insulator or a metal. The same ap- 
plies to similar temperature dependences observed, e.g., 
in Refs. [HI III III 



3. Weak-field regime 

As a matter of fact, the weak- field problem — whether 
or not there is an indefinite rise of the phase boundary 
as B — > — is a problem of the existence of a metal- 
insulator transition at B = and T = 0. In dilute 2D 
electron systems with low enough disorder, the resistiv- 
ity, p, strongly drops with lowering temperature provid- 
ing an independent way of facing the issue. Given strong 
temperature dependences of p, those with dp/dT > 
(dp/dT < 0) can be associated with a metallic (insu- 
lating) phase 00110 0. If extrapolation of 
the temperature dependences of p to T = is valid, the 
curve with dp/dT = should correspond to the metal- 
insulator transition. The fact that this method and the 
one, based on a vanishing activation energy combined 
with a vanishing nonlinearity of current-voltage curves 
when extrapolated from the insulating phase, give equiv- 
alent results strongly supports the existence of a true 
metal-insulator transition in zero magnetic field jl3 (see 
section III Dp . As long as in more-disordered 2D carrier 
systems the metallic (dp/dT > 0) behavior is suppressed 
(see, e.g., Refs. [IS EE E3, El E3, H3, El) or dis- 
appears entirely, it is definitely incorrect to extrapolate 
those weak temperature dependences of p to T — with 
the aim to distinguish between insulator and metal. Once 
one of the two methods fails, it remains to be seen how to 
verify the conclusion as inferred from the other method. 
This makes uncertain the existence of a true B = metal- 
insulator transition in 2D electron systems with high dis- 
order. 



4- Phase boundary oscillations 

The next important point is the oscillating behav- 
ior of the phase boundary that restricts the insulating 
phase with a xy = (see, e.g., Fig. It is worth not- 
ing that the oscillations persist down to the magnetic 
fields corresponding to the fillings of higher Landau lev- 
els, as indicated also by magnetoresistance oscillations 
0, E3 • The oscillation period includes the follow- 
ing elements. With decreasing magnetic field the low- 
est extended states follow the Landau level, float up in 
energy relative to its center, and merge with extended 
states in the next quantum level. The last element was 
absent in the original considerations 0,0,13113 lead- 
ing to discrepancies between experiment and theory. Re- 
cently, theoretical efforts have been concentrated on mod- 
ifications of the global phase diagram for the quantum 
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FIG. 4: (a) A sketch of the modified global phase diagram 
for the quantum Hall effect as expected from the mean-field 
approximation neglecting the Zeeman energy. The dashed 
line corresponds to the collapse of the exchange-enhanced spin 
splitting. Adopted from Ref. |86j|. (b) Numerical results for 
the phase diag ram within a tight-binding model. Adopted 
from Ref. 1031. 



Hall effect to reach topological compatibility with the ob- 
served metal-insulator phase diagram. It has been pre- 
dicted that the spin-up and spin-down extended states 
in the Landau level should merge |86ll87| (see Fig.|3Ja)). 
However, in view of lowest extended states, the topol- 
ogy of the phase diagram changes for the lowest Landau 
level only; besides, they do not float up before merg- 
ing. It has been verified that shifts of the extended 
states from the Landau level centers that are caused by 
disorder-induced mi xing of the Landau levels are small 
if. Within tight-binding models, 
an indication had been first obtained that the extended 
states disappear 0,0iE3i which caused some criticism 
toward the relevance of such a lattice model to the con- 
tinuum system [9?l 0|. After that, floating up of the 
extend ed states without merging has been foun d in stud- 
ies 0EO3IIO3- In contrast, Sheng and Weng pMll03| 
have obtained a merging of the extended states although 
without an oscillating behavior of the lowest extended 
states (see Fig.0Jb)). As for now, the effect of the phase 
boundary oscillations is still far from being described the- 
oretically. 
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Concluding this section, I would like to make a cou- 
ple of remarks on alternative ways for determining the 
metal-insulator phase boundary. An attempt was made 
to spot the phase boundary in the limit of B = using a 
criterion of a xy = e 2 /2h |l04j . However, that particular 
value of a xy has no special meaning as B — > 0. An idea 
was expressed to relate the minimum in the i nvers e com- 
pressibility to the metal- insulator transition 105]. How- 
ever, it has been recently shown that in zero magnetic 
field, this minimum lies at carrier densities well above 
the critica l den sity for the percolation metal-insulator 
transition |XQ6| - Particularly, in highly-disordered 2D 
carrier systems, its position may be close to that of the 
cross ing point of the resistivity at different temperatures 
|l05j | which formally yields overestimated densities for 
the metal-insulator transition because of suppression of 
the metallic behavior (see section Hi Dl . 



B. Similarity of the insulating phase and quantum 
Hall phases 

1. Method for comparison and consequences 
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About a decade ago, main attention was paid to the 
insulating phase at low electron densities as a possible 
candidate for the Wigner crystal. It was argued that 
its aforementioned reentrant behavior is a consequence 
of the competition between the quantum Hall effect and 
the pinned Wigner crystal [S |3jj . Another certain ar- 
gument was strongly nonlinear current-voltage charac- 
teristics in the insulating phase which were attributed 
to depinning of the Wigner crystal [40l Elj . Similar 
features of the insulating pha se in a 2D electron (near 
v = 1/5) EH ElHJpUp El El El and 2D hole 
(near v = 1/3) [^3, DDL uM systems in GaAs/AlGaAs 
heterostructures with relatively low disorder were also 
attributed to a pinned Wigner crystal which is inter- 
rupted by the fractional quantum Hall state. An al- 
ternative scenario was discuss ed in term s of percolation 
metal-insulator transition |57t Il07l Il08j . To distinguish 
between the two scenarios, the behavior of activation en- 
ergy and current-voltage characteristics in the insulating 
phase was stud ied a nd compared to that in quantum Hall 
phases (llUfillToflj. 

In contrast to the low-density insulating phase, the 
way of determining the current-voltage characteristics of 
the quantum Hall phases is different for Corbino and 
Hall bar geometries. In the former the dissipationless 
Hall current does not contribute to the dissipative cur- 
rent that is proportional to a xx , allowing straightforward 
measurements of current-voltage curves for all insulat- 
ing phases. In the latter the two current channels are 
connected through edge channels (see section HTC^ . and 
current-voltage characteristics correspond to quantum- 
Hall-effect breakdown curves. The dissipative backscat- 
tering current, /, that flows between opposite edge chan- 
nels is balanced by the Hall current in the filled Landau 



FIG. 5: Current-voltage characteristics in a low-disordered 
silicon MOSFET in B = 12 T at T ft! 25 mK for the low- 
density insulating phase at n s — 1.74 x 10 11 cm -2 (a) and the 
insulating phase with a xy h/e 2 = 1 at n s = 2.83 x 10 11 cm~ 2 
(b). In (b) the measured breakdown dependence Vxx(Isd) is 
converted into current- voltage characteristics (inset). From 
Ref. 0. 



levels associated with the longitudinal voltage, V xx . As 
long as a xx <C cr xy , the quantized value of a xy is a factor 
that allows determination of I = a xy V xx and the Hall 
voltage, V = Isd/&xy, from the experimental breakdown 
dependence of V xx on source-drain current, I s d- The de- 
pendence V(I) is a current- voltage characteristic, which 
is equivalent to the case of Corbino geometry [l£| (see 
Fig. E). Not only are the current-voltage curves simi- 
lar for all insulating phases, but they behave identically 
near the metal-insulator phase boundaries (see Fig.UJa)). 
The dependence of the critical voltage, V c , on the dis- 
tance from the phase boundary is close to a parabolic 
law 0, [57]]. The phase boundary position determined 
by a vanishing V c is practically coincident with that de- 
termined by a vanishing activation energy, E a , of elec- 
trons from the Fermi level Ep to the mobility edge, E c 
(see Fig. E^b)). The value E a is determined from the 
temperature dependence of the conduction in the linear 
interval of current-voltage curv es, which is activated at 
not too low temperatures |110| : note that it transforms 
into variable range hopping as T — ► (see below). The 
activation energy changes linearly with the distance from 
the phase boundary reflecting constancy of the thermo- 
dynamic density of states near the transition point (see 
also section lHDf) . The threshold behavior of the current- 
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FIG. 6: (a) Square root of the critical voltage as a function 
of electron density at the phase boundaries corresponding to 
a xy h/e 2 — 0, 1, and 2 in B = 12 T for a low-disordered 
2D electron system in silicon MOSFETs. (b) Behavior of 
the critical voltage and the activation energy near the phase 
boundary in B = 16 T. From Ref. [T§. 




B(T) 

FIG. 7: Metal-insulator phase diagram in a relatively low- 
disordered 2D electron system in GaAs/AlGaAs heterostruc- 
tures. The points corresponding to phase boundaries are ob- 
tained on two samples using a criterion of vanishing activation 
energy and vanishing nonlinearity of current-voltage curves as 
extrapolated from the insulati ng ph ase (circles and squares) 
and on one sample from Ref. |l!3t| using a cutoff criterion, 
°~xx = 10 MOhm, at a temperature » 25 mK as follows from 
the former method (diamonds). The solid lines are guides 
to the eye. Digits indicate a xy h/e 2 for different insulating 
phases. From Ref. |l6|. 



voltage characteristics is caused by the breakdown in the 
insulating phases. The breakdown occurs when the lo- 
calized electrons at the Fermi level gain enough energy 
to reach the mobility edge in an electric field, Vr/d, oy er 
a distance given by the localization length, L [la . Illl| : 

eV c L/d= \E C -E F \, (1) 

where d is the corresponding sample dimension. The 
values E a and V c are related through the localization 
length which is temperature independent and diverges 
near the transition as L(Ep) cx \E C — Ep\~ s with expo- 
nent s close to unity, in agreement with the theor etical 
value s = 4/3 in classical percolation problem |l!2j . The 
value of the localization length is practically the same 
near all metal-insulator phase boundaries, which indi- 
cates that even quantitatively, all insulating phases are 
very similar. Note that since the localization length in 
Eq. is small compared to the sample sizes, the phase 
boundary position determined by the diverging localiza- 
tion length refers to an infinite 2D system. As inferred 
from the vanishing of both E a and V c at the same point 
(see Fig. El^b)), possible shifts of the mobility threshold 
due to finite sample dimensions are small, which justifies 
extrapolations to the limit of L — > oo. 

The consequences of the method include (i) as long 



as there occurs no dramatic change in transport prop- 
erties, this excludes the pinned Wigner solid as the ori- 
gin for the insulating phase at low electron densities in 
available samples of low-disordered silicon MOSFETs; 
(ii) the metal- insulator phase diagram of Fig. [3Ja) is ver- 
ified and substantiated; (iii) the existence of a metal- 
insulator transition in zero magnetic field is supported 
(see section lll D(l : and (iv) the bandwidth of the extended 
states in the Landau levels is finite. All of these are 
also valid for relatively low-disordered 2D carrier sys- 
tems in GaAs/AlGaAs heterostructures with a distinc- 
tion that fractional quantum Hall phases are involved. 
Yet, the topology of the phase diagram remains un- 
changed including the oscillating behavior of the phase 
boundary that restricts the low-density insulating phase 
(see Fig. EJ). Additional confirmation of the percola- 
tion transition to the low-density insulating phase in 
GaAs/AlGaAs heterostructures was ob taine d by stud- 
ies of the high-frequency conductivity [ll3J_and time- 
resolved photoluminescence of 2D electrons |l!4j , as dis- 
cussed in Ref. 
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FIG. 8: Temperature dependence of the p xx peak width (n* — 
n s ) at half of the peak height counted from n* corresponding 
to v* = 2.5 (circles) and the behavior of the activation energy 
(triangles) in a low-disordered silicon MOSFET in B = 14 T. 
From Ref. 

2. Finite bandwidth of extended states 

It was predicted two decades ago that the localization 
length diverges as a power law at a sin gle energy, E* , 
which is the center of the Landau level |l!5l llltl Ill7j : 
L(E) oc \E — E*\~ s . An idea to check this predi ction 
based on low-temperature measurements of a xx |l!8j 
was qui ckly developed to a concept of single-parameter 
scaling |119( . It was suggested that the magnetoresis- 
tance tensor components are functions of a single vari- 
able that is determined by the ratio of the dephas- 
ing length, Ld{T) oc T~ p / 2 (where p is the inelastic- 
scattering-time exponent), and the localization length. 
The concept was claimed to be confirmed by measure- 
ments of temperature dependences of the peak width, 
AB, in p xx (or a xx ) and the maximum of dp xy /dB in 
a highly-disordered 2D electron system in InGaAs/InP 
heterostructures, which yielded AB oc T K , where k = 
p/2s « 0.4 |120| . Later, both deviations in the power law 
and different exponents in the range between k = 0.15 
and k = 1 were observed for other 2D carrier systems, 
different La ndau levels, and di f ferent disorder strengths 
|65j l68[ |6 9[ [l2ll Il22l H2I IT2II IT25I Il27l fl28l ll2flL 
113( 1 11311 11321 113.1 1134 113.1 I136L 11371 113ft Il3fll Il40j . 
Importantly, the scaling analysis of experimental data in 
question is based on two unverified assumptions: (i) zero 
bandwidth of the extended states in the Landau levels; 
and (ii) constancy of the thermodynamic density of states 
in the scaling range. If either assumption is not valid, this 
may lead at least to underestimating the experimental 
value of exponent k. 

The method of vanishing activation energy and van- 
ishing nonlinearity of current-voltage characteristics as 
extrapolated from the insulating phase shows that the 
former assumption is not justified. Moreover, measure- 
ments of the peak width in p xx as a function of temper- 
ature in low-disordered silicon MOSFETs yield a linear 
dependence which extrapolates as T — > to a finite peak 
width, accordingly [l^ (see Fig. |8J. Very similar tem- 
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FIG. 9: (a) Temperature dependence of the p xx peak width 
as determined by the maximum dp xy /dB at v* — 1.5 in a 
highly-disordered 2D electron system in GaAs/AlGaAs het- 
erostructures. Different symbols correspond to different runs. 
The dash ed lines are linear fits to the data. Adopted from 
Ref. |l2fil ]. (b) Localization length, determined by the high- 
frequency conductivity in the variable-range-hopping regime 
in a highly-disordered 2D electron system in GaAs/AlGaAs 
heterostructures, as a function of the filling factor deviation 
from v* = 2.5 (squares) and v* = 3.5 (circles) towards the 
v = 3 plateau. The dashed line is a linear fit that is expected 
from classical perc olatio n approach with exponent s = 4/3. 
Adopted from Ref. [lS^ . 



perature (and frequency) dependences were observed in 
highly-disordered 2D carr ier systems in GaAs/AlGaAs 
het erost ructures |l4llll42j| and Ge/SiGe heterostructures 
[6(1 Il43| . It is noteworthy that a similar behavior is re- 
vealed if the data from the publications, which claim the 
observation of scaling, is plotted on a linear rather than 
logarithmic scale (see, e.g., Fig. |§J; finite values of the 
peak width as T — > a r e eve n more conspicuous for the 
data of Refs. The reason for the ambi- 

guity is quite simple: within experimental uncertainty, 
it is difficult (on a logarithmic scale it is especially diffi- 
cult) to distinguish between sublinear/super linear fits to 
the data and linear fits which do not have to run to the 
origin. 

Although lack of data in most of the above experimen- 
tal papers does not allow one to verify the validity of 
both assumptions, it is very likely that there is no quali- 
tative difference between all of the discussed results. As 
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a matter of fact, they can be described by a linear (or 
weakly sublinear) temperature dependence with a finite 
offset at T = 0. Here comes an alternative and simple 
explanation of the temperature dependence of the peak 
width in p xx in terms of thermal broadening. Within 
percolation picture, if the activation energy E a ~ fesT, 
the conduction is order of the maximum a xx so that the 
value of ~ fcgT gives a thermal shift of the effective 
mobility edge corresponding to the a xx peak width [l5| . 
Despite the concept of thermal broadening has been ba- 
sically ignored in the literature in search for less trivial 
data interpretations, it looks as if no experimental results 
go beyond this favoring the concept of single-parameter 
scaling. Once the behavior of the localization length is 
not reflected by the temperature-dependent peak width 
in Pxx, no experimental support is provided for numeric 
calculations of the localization length which give a some- 
what larger exponent s « 2 compared to s = 4/3 in 
classical percolation problem (see, e.g., Ref. [l44] ) . The 
latter value of s as well as the behavior of the localiza- 
tion length in Fig. have been recently confirmed by 
measurements of the high- freq uency conductivity in the 
variable-range-hopping regime |138| (see Fig.^Jb)). 

Thus, the finding of finite bandwidth of the extended 
states in the Landau levels [THL IT^ |. which is in ob- 
vious contradiction to scaling arguments, is a firmly- 
established experimental result. Curiously, it has got no 
theoretical considerations for ten years. 

3. Hall insulator 

Deep in the insulating phases and at low temperatures 
the variable-range-hopping regime occurs in which the 
cond uctivity o xx is small compared to its peak value 
|l!2j | . In this regime it was predicted that the deviation, 
Act^ , of a xy from its quantized value in strong magnetic 
fields is much smaller than a xx oc exp(— (Tq/T) 1 / 2 ) |145| : 
A<T xy tx <jj x with exponent 7 « 1.5; note that this is in 
contrast to a straightforward linear relation in the acti- 
vation regime, as inferred from approximately the same 
behavior with temperature of the p xx peak width and the 
maximum of dp xy /dB. Later, finite p xy contrasted by di- 
verging p xx was found in calculations of the T = mag- 
netotransport coefficie nts in th e insulating phase with 
vanishing a xx and a xy |l46l Il47j . Such a behavior of p xx 
and p xy indicates a special quadratic relation between 
conductivities 

a xy oc a 2 xx . (2) 

Moreover, it wa s sho wn that p xy is close to the classical 
value (B/n s ec) |l48j providing arguments for the exis- 
tence of a Hall insulator phase |5Jl . 

Values p xy close to B/n s ec were experimentally found 
in the low-density insulating phase, deviations from the 
classical Hall line b eing attributed to po ssible admixture 
oi pxx @ El EH EHE EM EE1 EM EEi- Thus, the 
distinction of the Hall insulator phase from the quantum 
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FIG. 10: The induced voltage in a Corbino sample of a 
GaAs/AlGaAs heterostructure in up- and down-sweeps of the 
magnetic field. Also shown by straight lines are th e expected 
slopes for v = 2/3, 1, 2, 3, and 4. From Ref. |Tfi3 |. 

Hall phases — the absence of extended states below the 
Fermi level — becomes evident when expressed in terms 
of p xx and p xy . 

It was empirically found in low-disordered silicon 
MOSFETs that the lowest-filling-factor peak in a xx plot- 
ted in the (a^ ,,, o>y) plane is close to a semicircle centered 
at (e 2 /2/i,0) IM fiol fl54] . The semicircle law for the 
lowest-^ peak was reproduced in a highly-disordered 2D 
hole system in a Ge/SiGe quantum well |72i l73j . There, 
it was shown that the semicircle relation originates di- 
rectly from conductivity/resistivity tensor inversion 

/ e 2 \ 2 / e 2 \ 2 I- p , e 2 /h 

^ + {^-2h) ={2h) + pj x 7pl y ' (3) 

because the (narrow) a xx peak in question is located at 
filling factor just below v = 1 (see, e.g., Fig-EJ) where p xy 
is still close t o h/e 2 . Although this f inding is consistent 
with theories 155. ll5filT57i Il58l Fl59] . the semicircle law 
does not seem universal if highe r-^ peaks in a xx with 
different heights are involved |3^, Il23tll54l Il60| . 

C. Edge channel effects, direct measurements of 
the quantized Hall conductivity 

In a magnetically quantized 2D electron system, the 
Landau levels bend up at the sample edges due to the 
confining potential, and edge channels are formed whe re 
these intersect the Fermi energy (see, e.g., Ref. jl6l| V 
There rises a natural question whether the current in the 
quantum Hall state flows in the bulk or at the edges of 
the sample. Although the Hall conductivity a xy was not 
directly measured in early experiments on the quantum 
Hall effect, it seemed obvious that this value corresponds 
to p xy , in agree ment with the concept of currents that 
flow in the bulk |l62j |; that stands to reason that finite 
a xy would give evidence f or the exi stence of extended 
states in the Landau levels |l55l 1 1 6 1| . This concept was 
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challenged by the edge current model |l63j . In the lat- 
ter approach extended states in the bulk are not crucial 
and the problem of current distributions in the quantum 
Hall effect is reduced to a one-dimensional task in terms 
of transmission and reflection coefficients as defined by 
the backscattering current at the Fermi level between the 
edges. Importantly, if the edge current contributes signif- 
icantly to the net current, conductivity/resistivity tensor 
inversion is not justified, because the conductivities a xx 
and a xy are related to the bulk of the 2D electron system. 

To verify whether or not the Hall conductivity is quan- 
tized, direct measurements of a xy were necessary exclud- 
ing a possible shunting effect of the edge curre nts. Being 
equivalent to Laughlin's gedanken experiment |l64l Il65j |. 
such measurements were realized using the Corbino ge- 
ometry which allo ws separation of the bulk contribution 
to the net current pfi L Il67 l \)M Il70l Il7l l Il7lll7^ . 
A Hall charge transfer below the Fermi level between the 
coasts of a Corbino sample is induced by magnetic field 
sweep through the generated azimuthal electric field. If 
&xx 0, no discharge occurs allowing determination of 
the transferred charge, 

Q = a^TrrlffC^SB, (4) 

where r c g is the effective radius. The induced voltage, 
V = Q/C, which is restricted due to a large shunting ca- 
pacitance, C, changes linearly with magnetic field with a 
slope determined by a xy in the quantum Hall states until 
the dissipationless quantum Hall state breaks down (see 
Fig. 1101) . The fact that the quantization accuracy of a xy 
(about 1%) is worse compared to p xy may be attributed 
to non-constancy of the effective area in not very homo- 
geneous samples. Thus, the Hall current in the quantum 
Hall effect flows not only at the edges but also in the bulk 
of the 2D electron system through the extended states in 
the filled Landau levels. 

Apparently, the dissipative backscattering current in 
Hall bar samples should be balanced by Hall current in 
the filled Land au levels resulting in a longitudinal poten- 
tial drop |l74j . This point makes significant contribution 
to the edge current model. 

From an experimental viewpoint, all edge channel ef- 
fects proceed from slow equilibration (at macroscopic dis- 
tances) between the electrochemical potentials of differ- 
ent edge states including the state in the bulk. As long 
as such an equilibration occurs at the edges at the Fermi 
level, the applicability of the edge state model is justified. 
The approach accounts for effects observed in conven- 
tional transport experiments, which include the nonlocal 
resistance and the effects related to contacts/reservoirs 
(see, e.g., Ref. |175| ). However, particular potential pro- 
files at the edge and current distributions can be probed 
only using non-destructive spatially-re s olved imaging 
tech niques fl7fil Il77l Kj^L fl79L IjM fl8ll pi EH fl84 
IT MITMli!7Ul88Ul89Ul9rll^ note that 

many of the so-revealed inhomogeneous samples show 
quite good magnetotransport characteristics. Contrary 
to standard considerations of the edge channels in terms 



of skipping orbits for a confining potential that is sharp 
on the magnetic length, Ib = (hc/eB) 1 / 2 , it turns out 
that in most samples the potential profile at the edge is 
smooth and spans over much larger distances than lg. 
Edge regions corresponding to the scale of confining po- 
tential 10 /im were visualized in Hall photovoltage op- 
tical imaging experiments on standard Hall bar samples 
jl83j (see Fig. lTlT a'1'1. Since the Hall electric field is nearly 
consta nt, e ven if some field enhancement occurs near the 
edges |176| . the edge current contribution can be appre- 
ciable depending on a particular sample. 

For a soft confining potential, edge channels are also re- 
ferred to as compressible and incompressible strips whose 
spac ing is determined by the electron density gradient 
|195| . This is very similar to the long-standing phe- 
nomenon of Hall current pinch: given electron density 
gradients, the Hall current basically flows in narrow chan- 
nels (or incompressible strips) determined by the mini- 
mum a xx , their position in the sample being cont r olled , 
e.g., by magnetic field EH EH EH EH EH EH- 
When located at the edge, the pinch of Hall current be- 
comes identical with the subject of Ref. ,195]. The Hall 
current channels at the edge were imaged using a single- 
electron transistor as a probe for the local a xx |l84j (see 
Fig- ITTT b)). Applying a negative voltage to a side gate 
leads to a shift of the edge of the 2D electron system to- 
wards the probe, thereby creating a line scan of the local 
a xx across the sample edge. Vanishing a xx is indicated 
by enhanced fluctuations in the feedback signal, the feed- 
back circuit being used to keep the current through the 
single-electron transistor constant by controlling its volt- 
age relative to the 2D electron system. 

With respect to the preceding subsections, insignifi- 
cance of edge channel effects in transport experiments 
is verified in a usual way by coincidence of the results 
obtained in Hall bar and Corbino geometries. 



D. True zero-field metal-insulator transition, phase 
boundary in parallel magnetic fields 

As has been discussed above, the existence of extended 
states in quantizing magnetic fields is established by two 
independent experimental methods: (i) quantization of 
a xy ; and (ii) vanishing activation energy and vanishing 
nonlinearity of current-voltage characteristics as extrap- 
olated from the insulating phase. Theory is generally 
in agreement with this even though there are unresolved 
problems with finite bandwidth of the extended states 
in the Landau levels. In contrast, no extended states 
are expected in zero magnetic field, at least, for weakly- 
interacting 2D electron systems. The second experimen- 
tal criterion, however, results in an opposite conclusion 
although it does not have absolute credibility alone. To 
sort it out, further support by independent experimental 
verifications is needed. 

Another criterion is based on analysis of the temper- 
ature dependences of the resistivity in B = 0. Pro- 
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FIG. 11: (a) Image of the Hall photovoltage in a near surface GaAs/AlGaAs hetero junction at v — 2 (top) and a line scan 
taken horizontally through t he im age (bottom). The spot size is 5 pm across. The physical edges of the sample are indicated 
by vertical lines. From Ref. |l83l |. (b) Traces of the fluctuation part of the feedback signal vs the side-gate voltage for shifting 
the edge in a GaAs/AlGaAs heterostructure at different magnetic fields. Enhanced fluctuations in the feedback signal indicate 
vanishing a xx , which corresponds to Hall current channels (incompressible strips) for v = 1, 2, 3, and 4. Adopted from 
Ref. [HI. 



vided these are strong, those with positive (negative) 
derivative dp/dT are indicative of a metal (insulator) 

EmiHlzlIzi; note that in the vicinity of the tran- 
sition, p(T) dependences obey the scaling law with expo- 
nent k « 1, which is consistent with the concept of ther- 
mal broadening/shift by the value ~ ksT of the effective 
mobility edge in the insulating phase (see section lTlBfl . If 
extrapolation of p(T) to T = is valid, the critical point 
for the metal- insulator transition is given by dp/ dT = 0. 
In a low-disordered 2D electron system in silicon MOS- 
FETs, the resistivity at a certain electron density shows 
virtually no te mperature dep endence over a wide range of 
temperatures |2Q3 EMI EE>3 (see Fig. [Ha)). This curve 
separates those with positive and negative dp /dT nearly 
symmetrically at temperatures above 0.2 K 203]. As- 
suming that it remains flat down to T = Q, one obtains 
the critical point n c which corresponds to a resistivity 
p« 3/i/e 2 M 

Recently, these two criteria have been applied simul- 
taneously to the 2D metal-insulator transition in low- 



disordered silicon MOSFETs [73, l202j . In zero magnetic 
field, both methods yield the same critical density (see 
Figs. I12f b) and I13f b)). Since one of the methods is 
temperature independent, this equivalence strongly sup- 
ports the existence of a true metal-insulator transition 
in 5 = 0. This also adds confidence that the curve 
with zero derivative dp/dT will remain flat (or at least 
will retain finite resistivity value) down to zero temper- 
ature. Additional confirmation in favor of true zero-field 
metal-insulator transition is provided by magnetic mea- 
surements as described in the next section. 

In the presence of a parallel magnetic field, B» , the out- 
come is very different. With increasing parallel field the 
transition point, n c (Bu), determined from the vanishing 
nonlinearity and activation energy, shifts approximately 
linearly to higher electron densities, saturating above a 
critical field, B c w 3 T, at a constant value which is ap- 
proximately 1.5 times higher than that in zero field. Note 
that a similar suppression of the metallic state was ob- 
served using a cutoff criterion, p = 100 kOhm |l08j |. In 
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FIG. 12: (a) Resistivity as a function of temperature at 
different electron densities in a low-disordered silicon MOS- 
FET. The inset sho ws th e middle curve on an expanded scale. 
Adopted from Ref. |201| . (b, c) Temperature dependence of 
the resistivity of a low-disordered silicon MOSFET at differ- 
ent electron densities near the metal-insulator transition in 
zero magnetic field (b) and in a parallel magnetic field of 4 T 
(c). The electron densities are indicated in units of 10 11 cm -2 . 
From Ref. El. 



the metallic phase the saturation of the resistance with 
parallel field signals the onset of full spin polarization 
of the 2D electrons, as inferred from an analysis of the 
Shubnikov-de Haas oscillations in tilted magnetic fields 
|204L 12051 l206| . One expects that the 2D electron system 
is spin polarized at parallel fields B\\ > B c , and that the 
observed phase boundary shift is a spin effect. At the so- 
determined critical density n c (B\\), the exponential diver- 
gence of the resistivity as T — > ends, although dp/dT 
remains negative at least for Bn > B c (see Figs. I12f c) 
and E{ b)). In sharp contrast with the zero-field case, 
not only are the p(T) curves in the field nonsymmetric 
about the middle curve in Fig. I12f . but all of them 
have negative derivatives dp/dT in the entire tempera- 
ture range, although the values of p are comparable to 
those in the B = case. The metallic (dp/dT > 0) tem- 
perature dependence of the resistivity, observed at higher 
electron densities in parallel magnetic fields, is weak so 
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FIG. 13: (a) Current-voltage characteristics in zero mag- 
netic field at ~ 30 (solid line) and 211 mK (dashed line) for 
the same silicon MOSFET as in Fig. EH b, c); note that the 
threshold voltage is practically independent of temperature. 
An Arrhenius plot of the resistivity in the insulating phase 
is displayed in the inset for different values of Bu and n s . 
From Ref. [77|. (b) Activation energy and square root of the 
threshold voltage as a function of electron density in zero 
magnetic field (circles) and in a parallel magnetic field of 4 T 
(diamonds). The critical densities correspond to the dashed 
lines in Fig. I12f b. c). From Ref. (c) Resistivity ver- 

sus temperature in a strongly-disordered silicon MOSFET at 
the following electron densities: 3.85, 4.13, 4.83, 5.53, 6.23, 
7.63, 9.03, 10.4, 11.8, 13.2, 16.0, 18.8, 21.6, 24.4, 30.0, and 
37.0 x 10 11 cm -2 . The p(n s ) isotherms are shown in the in- 
set. Adopted from Ref. 



that the derivative method does not yield a critical den- 
sity at -B|| > B c . Its failure leaves uncertain the existence 
of a true metal-insulator transition in a parallel magnetic 
field. 

A very similar conclusion holds for 2D electron sys- 
tems with higher disorder in zero magnetic field (see sec- 
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FIG. 14: Shubnikov-de Haas oscillations in a low-disordered silicon MOSFET at a temperature ~ 40 niK for high (a) and 
low (b) electron densities, (c, d) Evolution of the Shubnikov-de Haas oscillations with electron density in two ranges of filling 
factors. The curves are arbitrarily shifted along the t/-axis for clarity. From Ref. [2^ |. 



tion 111 Aj l . In this case the metallic (dp/dT > 0) behav- 
ior is also suppressed [IE E3, E3, Ell 113, E3, E3, E3| or 

disappears entirely, and extrapolation of the weak p(T) 
dependences to T = is not justified invalidating the 
derivative criterion for the critical point for the metal- 
insulator transition (see Fig. I13f c)). This is notewor- 
thy that owing to its simplicity, the derivative method 
is widely used for describing metallic (dp/dT > 0) and 
insulating (dp/dT < 0) temperature dependences of re- 
sistance in a restricted temperature range. To avoid con- 
fusion with metallic and insulating phases, however, one 
should employ alternative methods for determining the 
metal-insulator transition point. Such methods including 
a vanishing activation energy and noise measurements 
have been appl i ed to highly-disordered 2D carrier sys- 
tems 202[ 12071 l208j . Being similar, they yield lower 
critical densities for the metal-insulator transition com- 
pared to those obtained using formally the derivative cri- 
terion. This simply reflects the fact that the metallic 
(dp/dT > 0) behavior is suppressed, the critical den- 
sity n c increasing naturally with disorder strength (see 
Fig.Q}. 



III. MANY-BODY PHENOMENA IN DILUTE 
2D ELECTRON SYSTEMS 



The resistivity drop with decreasing temperature in 
a low-disordered 2D electron system in silicon MOS- 
FETs in zero magnetic field being strong compared to the 
metallic p(T) expe cted from the temp erature-dependent 
screening theories |209t l210t l21ll |212| , it was attributed 
to a manifestation of strong electron-electron interactions 
|l9). Recently, the underlying physics of the effect has 
been clarified. At low electron densities, a strongly en- 
hanced ratio gm of the spin and the cyclotron splittings 
has been found indicating that the 2D system behaves 
well beyond the weakly interacting Fermi liquid . Ex- 
perimental results have also shown that it is the effec- 
tive mass that increases sharply at low electron densi- 
ties and is related to the anomalous rise of the resistiv- 
ity with temperature |28|. In view of T = quantum 
phase transitions, while for the insulating phase a transi- 
tion is signaled by the diverging localization length, the 
interaction-enhanced mass may be a similar indicator for 
the metallic phase. 
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A. Strong increase of the product gm near the 
metal-insulator transition, possible ferromagnetic 
transition 



1. Beating pattern of Shubnikov-de Haas oscillations 



Electron-electron interactions give rise to a renormal- 
ization of the Fermi-liquid parameters including the ef- 
fective mass and g factor Q . Tracing Shubnikov-de Haas 
oscillation minima in a 2D electron system in tilted mag- 
netic fields, it is easy to determine the ratio gm of spin 
and cyclotron splittings which is proportional to the spin 
susceptibility \- I n the range of high electron densities 
> 2 x 10 11 cm~ 2 in silicon MOSFETs, moderate en- 
hancements of gm by a factor of < 2.5 were observed 
[204L l213t l214j , which is consistent with the concept of 
weakly interacting Fermi liquid. 

At low electron densities in low-disordered sili- 
con MOSFETs in perpendicular magnetic fields, the 
Shubnikov-de Haas oscillation minima corresponding to 
the cyclotron splittings (y — 4, 8, 12, 16,...) were found 
to disappear as the electron density is reduced (see 
Fig. I14fl . Disregarding the minimum for the valley split- 
ting at v = 1, only minima corresponding to the spin 
splittings (y = 2, 6, 10, 14,...) remain close to the metal- 
insulator transition which occurs in the studied samples 
at n c w 8 x 10 10 cm -2 . These results show that as one 
approaches the metal-insulator transition, the cyclotron 
gaps (which are equal to the difference between the cy- 
clotron and spin splittings, ignoring the valley splitting) 
become smaller than the spin gaps and eventually van- 
ish. The condition for vanishing is coincidence of the spin 
and the cyclotron splittings, or gm/2m e = 1 (where m e 
is the free electron mass), which is higher by more than 
a factor of 5 than the value of this ratio in bulk silicon, 
gmj1m e = 0.19. The effect cannot be explained in the 
framework of the many-body en hancement of spin gaps 
in a perpendicular magnetic field pl^ El^ 121^ 121^ l219j 
because the disappearance of the cyclotron gaps in a wide 
range of magnetic fields would require an enhanced g fac- 
tor which is independent of magnetic field. This implies 
that the product gm is nearly field-independent and ap- 
proximately equal to its many-body enhanced zero-field 
value (see section IlII Bfl . Thus, the spin susceptibility 
X oc gm is strongly enhanced near the metal-insulator 
transition. 

Similar experiments in tilted magnetic fields cannot 
provide high accuracy in determining the behavior of 
the renormalized gm at low electron densities because 
there are too few Shubnikov-de Haas oscillations near the 
metal-insulator transition. High accuracy was attained in 
experiments on the parallel- field magnetotransport. 



2. Scaling the parallel-field magnetoresistance and other 
methods 

As the thickness of the 2D electron system in silicon 
MOSFETs is small compared to the magnetic length 
in accessible fields, the parallel magnetic field couples 
largely to the electrons' spins while the orbital effects are 
suppressed. The resistance in dilute silicon MOSFETs 
was found to be isotropic with respect to in-plane mag- 
netic field and rise steeply with the field saturating to a 
constant value a bove a critical fi eld B c which depends on 
electron density [220l l22ll l222j (see Fig. E^a)). As has 
been mentioned in section III Dl the saturation field B c 
corresponds to t he onset of fu ll spin polarization of the 
electron system po l l205l l206j . 

In a low-disordered 2D electron system in silicon MOS- 
FETs, it was found that the normalized magnetoresis- 
tance, measured at different electron densities in the 
low-temperature limit in which p(B^) becomes temper- 
ature independent, collapses onto a single curve when 
plotted as a function of B\\/B c , where the scaling pa- 
rameter B c is normalized to correspond to the satura- 
tion/polarization field [23| (see Fig. lT5T b. c)). The scaling 
breaks down when one approaches the metal-insulator 
transition where the magnetoresistance depends strongly 
on temperature even at the lowest experimentally achiev- 
able temperatures. Note that the observed scaling de- 
pendence is described reasonably well by the theoreti- 
cal dependence of p/p(0) on the degree of spin polar- 
ization, £ = gmpBB\\/ Trh 2 n s = B\\/B c , based on the 
spin-polarization-dependent screening of a random po- 
tential |223j| . The field B c is proportional, with high 
precision, to the deviation of the electron density from 
its critical value: B c oc (n s — n c ) (see Fig. llCf a)). The 
procedure used provides high accuracy for determining 
the functional form of B c (n s ), even though the abso- 
lute value of B c is determined not so accurately. Since 
the strong increase of the product gm at low electron 
densities (see Fig. I16f b)). which follows from the B c (n s ) 
dependence, is in agreement with the enhanced gm ob- 
tained by Shubnikov-de Haas oscillations, the band tail 
of localized electron states is small and the clean limit 
occurs. Therefore, the tendency for B c to vanish at a 
finite electron density n x close to n c gives evidence in fa- 
vor of the existence of a ferromagnetic transition in this 
electron system indicating that the metal-insulator tran- 
sition is driven by interactions p3 |. It signifies that on 
the phase diagram of Fig.^ the vicinity of either tricrit- 
ical point is reached (see section UlI Cll . 

A similar conclusion about possible spontaneous spin 
polarization was drawn based on a scaling of magneto- 
conductivity data in similar samples at different electron 
densities and temperatures j24[. Results for the strongly 
enhanced gm were corroborated in detailed studies of 
Shubnikov-de Haas oscillations in dilute silicon MOS- 
FETs with higher disorder in tilted magnetic fields 224] 
(see Fig. I16f c)). The agreement between all three sets 
of data is remarkable, especially if one takes into ac- 
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FIG. 15: (a) Resistivity versus parallel magnetic field mea- 
sured on a low-disordered silicon MOSFET. Different symbols 
correspond to gate voltages from 1.55 to 2.6 V, or to densi- 
ties from 1.01 to 2.17 x 10 11 cm" 2 . Adopted from Ref. |22ljl . 
(b) Low-temperature magnetoresistance of a low-disordered 
2D electron system in silicon MOSFETs in parallel magnetic 
fields at different densities above n c for the B = metal- 
insulator transition. From Ref. p^. (c) Scaled curves of the 
normalized magnetoresistance vs B^/B c . The electron den- 
sities are indicated in units of 10 11 cm -2 . Also shown by a 
dashed line is the nor malize d magnetoresistance calculated by 
Dolgopolov and Gold ^22^. From Ref. [23| . 



count that different groups used different methods, differ- 
ent samples, and different field/spin-polarization ranges 
|2fij . This also indicates that the electron density n x ~ 
8 x 10 10 cm~ 2 is sample independent, in contrast to the 
critical density n c for the metal-insulator transition. Ob- 
viously, for the spin susceptibility x °t 9 m to diverge at 
n s = n x , the extrapolation of B c (n s ) to zero must be 
valid. To verify its validity, accurate data at lower den- 
sities, lower temperatur es, and o n much less disordered 
samples are needed pl Ezfl l226j . 

Thermodynamic investigations of the spin suscep- 
tibility, based on measurements of the chemical po- 
tential change with parallel magnetic field, dfi/dBu, 
were performed in highly-disordered silicon MOSFETs, 
as inferred from the consider ably higher densities for 
the metal-insulator transition |227| (see Fig. I16f d). cf. 
Fig. Ilfif c)). As compared to the clean regime, the ob- 
tained dependence of the polarization field B c on n s 
in Fig. I16f d) is shifted to appreciably higher electron 
dens i ties caused b y local moments in the band tail 
[221 12271 122II22H I. The band tail effects thus become 
crucial in parallel-field experiments on highly-disordered 
2D electron systems. 



3. Other 2D carrier systems 

In dilute GaAs/AlGaAs heterostructures, a similar 
enhancement of the spin susceptibility at low electron 
densities was found by an analysis of the Shubnikov- 
de Haas oscillations [23 (see Fig. UTT a)). The thick- 
ness of the 2D carrier system in GaAs/AlGaAs het- 
erostructures is relatively large, which leads to an in- 
crea se of the effective mass w ith parallel magnetic field 
[H EU EH iai- As a result, the value of 
polarization field obtained from the parallel-field mag- 
netoresistance becomes strongly reduced as the electron 
density increases. Disregarding this reduction, both data 
sets determine the B c (n s ) dependence whose critical be- 
havior is not so evident, possibly because the lowest ex- 
perimentally reached densities are still too high. Due 
to the lower effective mass, higher dielectric constant, 
and the absence of valley degeneracy, the same inter- 
action strength r* = E ee /Ep in the 2D electron sys- 
tem in GaAs/AlGaAs heterostructures is expected to bc 
achieved at densities about two orders of magnitude lower 
than in silicon MOSFETs. Therefore, the critical re- 
gion is expected at densities n s < 10 9 cm -2 which have 
not yet been accessed in currently available samples of 
GaAs/AlGaAs heterostructures. 

The orbital effects in parallel magnetic fields can be 
avoided by using narrow quantum wells. A 2D electron 
system in narrow AlAs quantum wells is similar to that 
in silicon MOSFETs, excep t that in the former the valley 
degeneracy is absent |235j . The critical region expected 
at densities n s < 2 x 10 10 enr 2 is exceeded strongly by 
the lowest accessible electron densities in narrow AlAs 
quantum wells with high disorder (see Fig. 117( b)). Note 
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FIG. 16: (a) Dependence of the field B c on electron density in a low-disordered silicon MOSFET. The dashed line is a linear fit. 
From Ref. |23|. (b) The product gm versus electron density obtained from the data for B c . From Ref. [23|. (c) Polarization field 
as a function of electron density obtained by different groups 23, 24, 224]. The critical density n c in the samples of Ref. .23] is 
indicated. From Ref. U^l . (d) Polarization field vs electron density in a highly-disordered silicon MOSFET as determined from 
dfi/dB^ — at 0.2 K. The dashed line is a linear fit. The evaluated position of the metal-insulator transition is indicated. The 
derivative d/j,/dB^ as a function of el ectron density in a parallel magnetic field of 9 T at a temperature of 100 mK is displayed 
in the inset. Adopted from Ref. |227| . 



that the data points in the insulating ph ase reflect the 
phy sics of local moments in the band tail [222L 12271 12281 
[229J, which is different from that of the metallic phase. 



Being very similar to silicon MOSFETs, a 2D elec- 
tron system in Si/SiGe quantum wells is different by 
the higher dielectric constant and the presence of a 
spacer. Moreover, it is distinguished from other sys- 
tems by its remote-doping scattering as indi cated by the 
small parallel-field magnetoresistance |236j . A similar 
increase of the spin susceptibility at low densities was 
observed in this electron system, the lowest achievable 
densities also being well abo ve t he expected critical re- 
gion at n s < 4 x 10 10 cm~ 2 |237f . So, in all the studied 
dilute 2D electron systems other than silicon MOSFETs, 
too high disorder and, hence, too high densities for the 
metal-insulator transition (see Fig. ^) mask possible crit- 
ical behavior of the spin susceptibility. 



B. Determining separately the effective mass and g 

factor 

1. Slope of the metallic temperature dependence of 
conductivity in zero magnetic field 

The strong enhancement of the spin susceptibility x c< - 
gm at low electron densities can be caused in principle by 
an increase of either g or m or both. The effective mass 
and g factor were determined separately using the recent 
theory of temperature-dependent corrections to conduc- 
tivity due to electron-electron interactions |238| |. Note 
that its main advantage com pared to the tempe rature- 
dependent screening theories |209l 

IH3 HO, IH3 is that 
spin exchange effects arc treated carefully in the new the- 
ory. At intermediate temperatures, the predicted <r(T) is 
a linear function 

^l = l-A*k B T, A* = - ( 1 + 8 /^ , (5) 
do Trn z n s 
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FIG. 17: (a) Data for B c as a function of electron density in a 
dilute GaAs/AlGaAs heterostructure obtained by Shubnikov- 
de Haas oscillations (squares) and parallel-field magnetoresis- 
tance (triangles). The power-law fit B c oc n\ A (solid line) is 
compared to the linear fit (dashed line). Also shown is the 
evaluated position of the metal-insulator transition. Adopted 
from Ref. |27| . (b) Dependence of B c on electron density 
obtained by measurements of the parallel-field magnetoresis- 
tance in highly-disordered samples of narrow AlAs quantum 
wells. The solid curve represents the quantum Monte Carlo 
calculation for a disorder-free 2D electron system The 
dashed line is a linear fit. The evaluated position of th e me tal- 
insulator transition is indicated. Adopted from Ref. |235l] . 



where <7o is the value obtained by extrapolating the lin- 
ear interval of the a(T) dependence to T = and the 
factor of 8 is expected for temperatures lower than the 
valley splitting in silicon MOSFETs. The slope, A*, is 
determined by the Fermi liquid constants Fq and F* that 
define the renormalization of g and m: g/go — 1/(1 + Fq) 
and m/mb = 1 + F{ ■ Using these relations one obtains 
both g and m from the data for the slope A* and the 
product gm psf . 

For sufficiently small deviations |it/(Tq — 1|, the depen- 
dence of the normalized conductivity er/<7o on tempera- 
ture at different electron densities above the critical den- 
sity n c for the metal-insulator transition is linear over 
a wide enough interval of temperatures (see Fig. HHT a)). 
The inverse slope 1/A* and the value hbB c are close 
to each other in a wide range of electron densities (see 
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FIG. 18: (a) The temperature dependence of the normalized 
conductivity at different electron densities (indicated in units 
of 10 11 cm -2 ) in a low-disordered silicon MOSFET above 
the critical electron density for the metal-insulator transition. 
The dashed lines are fits to the linear interval of the depen- 
dence, (b) The inverse slope 1/A* (circles) and the polariza- 
tion field B c (diamonds) as a function of electron density. The 
dashed lines are linear fits which extrapolate to the critical 
density for the metal-insulator transition, (c) The effective 
mass and g factor versus electron density determined from 
an analysis of the temperature-dependent conductivity and 
parallel-field magnetoresistance. The dashed lines are guides 
to the eye. From Ref. H|. 
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Fig. HBTb)). Moreover, the low density data for 1/A* 
are approximated well by a linear dependence which ex- 
trapolates to the critical density n c in a way similar to 
the behavior of B c . This finding immediately points to 
approximate constancy of the g factor at low electron 
densities, according to the functional form of the slope 
A* in Eq. Q. 



Renormalizations g/go and mjra^ versus electron den- 
sity determined from this analysis corroborate earlier re- 
sults at high densities but are striking in the limit of low 
electron densities (see Fig. ll8f cl). In the high n s region, 
the enhancement of both g and m is relatively small, both 
values increasing slightly with decreasing electron density 
in agreement with earlier data 239]. Also, the renor- 
malization of the g factor is dominant compared to that 
of the effective mass, consistent with theoretical studies 
|24ft l24ll l242j . In contrast, in the low n s region, the 
renormalization of the effective mass increases sharply 
with decreasing density while the g factor remains nearly 
constant. Hence, it is the effective mass, rather than the 
g factor, that is responsible for the drastically enhanced 
spin susceptibility near the metal-insulator transition. 
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Normally, no difference is assumed between the inter- 
action parameter r* = E ee / 'Ep and the Wigner-Seitz ra- 
dius r s = 1/ (-Kris) 1 / 2 a,B- The finding of the strongly 
enhanced effective mass breaks the equivalence of r* and 
r s because they are connected through the n s dependent 
mass: r* = 2(m/mb)r s (where the factor of 2 comes from 
the valley degeneracy in silicon MOSFETs). Therefore, 
as one approaches the metal-insulator transition in low- 
disordered silicon MOSFETs, the interaction parameter 
r* grows much more rapidly than r s reaching the values 
r* s > 50 H. 

In similar experimental verifications of the theory 
238] on different 2D carrier systems, much higher val- 
ues Fq than the expected limit Fq — —1 for the 
Stone r ins t ability were fou nd for the metallic slope of 
cr(T) HI HI HH HH. The relatively small en- 
hancement of the g factor indicates that the spin ex- 
change effects are not very pronounced. From an ex- 
perimental point of view, this raises a problem of com- 
parison between different theories. Formally, even if 
there are possible uncertainties in the theoretical coef- 
ficients, both the temp erature-depende nt sc reening theo- 
ries [20a, EH and the theory 12381 des cribe rea- 
sonably well available experimental data |24il |248| . To 
discriminate between these two, more detailed compari- 
son with experiment is needed. Note that in all theories, 
the slope A* in Eq. © is proportional to the effective 
mass, so the conclusion about the strongly enhanced ef- 
fective mass at low electron densities in low-disordered 
silicon MOSFETs is basically independent of a particu- 
lar theory. 
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FIG. 19: (a) Change of the amplitude of the weak-field 
Shubnikov-de Haas oscillations in a low-disordered silicon 
MOSFET with temperature at n s = 1.17 x 10 11 cm~ 2 for 
oscillation numbers u = 10 (circles) and v = 14 (squares). 
The value of T for the v — 10 data is divided by the factor 
of 1.4. The solid line is a fit using Eq. ©. (b) Dependence 
of the effective mass on electron density determined from an 
analysis of Shubnikov-de Haas oscillations (circles) and from 
an analysis of o(T) and p(B«) (dashed line), (c) The effective 
mass versus the degree of spin polarization for the following 
electron densities in units of 10 11 cm -2 : 1.32 (circles), 1.47 
(squares), 2.07 (diamonds), and 2.67 (triangles). The dashed 
lines are guides to the eye. From Ref. [2S| . 
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2. Temperature dependent amplitude of the weak-field 
Shubnikov-de Haas oscillations 

The claim about the strong increase of the effective 
mass was verified based on analysis of the temperature 
dependence of the Shubnikov-de Haas oscillations. The 
method is similar to that used by Smith and Stiles |214j 
but is extended to much lower electron densities and tem- 
peratures [29|. In the low temperature limit, the p(T) 
dependence saturates, and the Lifshitz-Kosevich formula 
with constant Dingle temperature for the weak-field os- 
cillation amplitude for the normalized resistance 

A{T) 2TT 2 k B T/Hu) c 
A ~ sinh(27r 2 fc s T/fiij c )' 
A = 4exp(-2ir 2 k B T D /htu c ) (6) 

(where oj c = eB^/mc is the cyclotron frequency and 
Try is the Dingle temperature) describes damping of the 
Shubnikov-de Haas oscillations with temperature (see 
Fig. H9f a)). The effective mass as a function of electron 
density determined by this method agrees well with the 
data obtained by the procedure described in the preced- 
ing section (see Fig. Iiyf b)'). The agreement between the 
results obtained using two independent methods supports 
the validity of both and justifies applicability of Eq. © 
to the strongly interacting 2D electron system in silicon 
MOSFETs. 

To probe a possible contribution from the spin ex- 
change effects to the effective mass enhancement, a paral- 
lel magnetic field component was introduced to align the 
electrons' spins. Within the experimental accuracy, the 
effective mass does not depend on the degree of spin po- 
larization f = (B\ + B 2 f/ 2 /B c (see Fig.QjUc)). There- 
fore, the m(n s ) dependence is robust, and the origin of 
the mass enhancement has no relation to the electrons' 
spins and exchange effects. 

A similar analysis of the Shubnikov-de Haas oscilla- 
tions in dilute silicon MOSFETs at high temperatures 
T > 0.3 K, where the low-density resistivity (and, hence, 
Td) depends strongly on temperature, allows an evalua- 
tion of the effective mass as well as th e q fa ctor which is 
calculated from the known value of gm |224j (see Fig. 121)11 . 
The two data sets are obtained based on an assumption 
of a temperature-independent Dingle temperature and 
that of a Dingle temperature that increases linearly with 
temperature. Too large uncertainty of the results makes 
it impossible to establish which value (either g or m or 
both) is responsible for the strong enhancement of the 
spin susceptibility. Note that an attempt to improve the 
evaluation of the effective mass by justifying the applica- 
tion of the Lif shitz -Kosevich formula with temperature 
dependent Try [249j would lead on the contrary to bigger 
deviations of the estimated values of mass from the data 
obtained in the low temperature limit (cf. Figs. I19f 
and !2U|) . 

An analysis of the temperature dependent amplitude 
of Shubnikov-de Haas oscillations in a dilute 2D electron 
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FIG. 20: The effective mass and g factor in a dilute silicon 
MOSFET as a function of r s . Different symbols correspond 
to two different assumptions for evaluating m from an anal- 
ysis of the high-temperature Shubnikov-de Haas oscillations: 
a temperature-independent Td (open circles) and a Dingle 
temperature that increases linearly with temperature (solid 
circles). The solid a nd d ashed lines in (a) are polynomial fits. 
Adopted from Ref. ^H - 



system in narrow AlAs quantum wells yielded moderate 
enhancements of the effective mass as we ll as the g fac- 
tor determined from the known grn 235]. The observed 
behavior of g and m is similar to that found at high 
electron densities in silicon MOSFETs. This indicates 
that the valley origin of the strongly enhanced effective 
mass at low electron densities in silicon MOSFETs is not 
very likely, even though the lowest accessible densities 
in narrow AlAs quantum wells are still too high. Inter- 
estingly, the observed values of g/go in the limit of high 
electron densities in AlAs quantum wells exceed appre- 
ciably the value of g/go = 1 as well as those in silicon 
MOSFETs. The increase of the spin susceptibility with 
the strain-induced valley polarization, ob serve d at high 
electron densities in AlAs quantum wells |250j |. is likely 
to be connected with the increase in the g factor. 



3. Spin and cyclotron gaps in strong magnetic fields 

The results for the strong enhancement of the effec- 
tive mass are also consistent with the data for spin and 
cyclotron gaps obtained by magnetocapacitance spec- 
troscopy. The method is based on determination of the 
chemical potential jumps in a 2D electron system when 
the filling factor traverses the gaps in the spectrum. A 
dip in the magnetocapacitance at integer filling factor is 
directly related to a jump of the chemical potential across 
a corres ponding gap in the spectrum of the 2D electron 
system |25lll25^ : 
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FIG. 21: For silicon MOSFETs, difference of the normalized 
values of cyclotron and spin gaps in a perpendicular magnetic 
field versus electron density. The level width contribution is 
indicated by systematic error bars. Also shown for compar- 
ison is the value (m e /m — g) determined from the da ta of 
Ref. [H (dashed line), Ref. [2<| (dotted line), and Ref. plH| 
(dash-dotted line) as well as using the band electro n ma ss and 
the g factor in bulk silicon (solid line). From Ref. |254f . 



where Co is the geometric capacitance and A g is the sam- 
ple area. The chemical potential jump is determined by 
integrating the magnetocapacitance over the dip in the 
low temperature limit where the magnetocapacitanc e sat - 
urates and becomes independent of temperature |253| . 
Note that standard measurements of activation energy 
yield a mobility gap which may be different from the 
gap in the spectrum. This is a serious disadvantage as 
compared to the direct method of magnetocapacitance 
spectroscopy. 



The g factor determined by this method is close to its 
value i n bu lk silicon and does not change with the filling 
factor 254], in disag r eem e nt w i th the the ory of exchange- 
enhanced gaps pTa . 121^. 12171 l219j . The cyclotron 
splitting corresponds to the effective mass that is strongly 
enhanced at low electron densities (see Fig.l21fl. Thus, in 
strong magnetic fields, the spin exchange effects are still 
not pronounced. 



It is worth noting that in contrast to the g factor, 
the valley gap is strongly enhanced at the lowest fill- 
ing factors v = 1 and v = 3 and oscillates with v 
|253l |254[ . This is similar to the beha vior of th e spin 
gap in GaAs/AlGaAs heterostructures 255, 25G], both 
of the gaps increasing linearly with perpendicular mag- 
netic field. 



C. Wigner crystal or ferromagnetic Fermi liquid, 
theoretical approaches 

As has been mentioned above, the experimental re- 
sults obtained in low-disordered silicon MOSFETs indi- 
cate that on the metallic side the metal-insulator transi- 
tion is driven by interactions. In contrast, on the insulat- 
ing side this is still a classical percolation transition with 
no dramatic effects from interactions. One concludes that 
the vicinity of either tricritical point on the phase dia- 
gram in Fig. n is reached. This is consistent with the 
fact that the interaction parameter r* at low electron 
densities exceeds the theoretical estimate for the onset 
of Wigner crystallization, even though it is not yet clear 
whether or not electron crystallization expected in the 
low density limit is preceded by an intermediate phase 
like ferromagnetic Fermi liquid. 

To face the problem, two approaches have been for- 
mulated. The first one exploits the Fermi liquid model 
extending it to relatively large r*. Its outcome is that 
the renorma l izati on of g is large compared to that of 
m |24fl l24ll l242j . In the limiting case of high r * , one 
may expect a divergence of the g factor that corresponds 
to the Stoner instability. These predictions are in con- 
tradiction to the experimental data. Firstly, the dilute 
system behavior in the regime of the strongly enhanced 
susceptibility — close to the onset of spontaneous spin 
polarization and Wigner crystallization — is governed by 
the effective mass, rather than the g factor, through the 
interaction parameter r * . And secondly, the insensitivity 
of the effective mass enhancement to spin exchange ef- 
fects cannot be accounted for. This discrepancy reduces 
somewhat the chances for the occurrence of a ferromag- 
netic Fermi liquid that precedes electron crystallization. 
In principle, should the spin exchange be small, the spin 
effects may still come into play closer to the onset of 
Wigner crystallization where the Fermi energy may con- 
tinue dropping as caused by mass enhancement. 

The other theoretical approach is not based on Fermi 
liquid. In analogy with He 3 , it was predicted the exis- 
tence of an intermediate phase between the Fermi liquid 
and the Wigner crysta l, ca used by a partial separation 
of the uniform phases [257]. It was also predicted that 
near the crystallization point, the renormalization of 777 
is dominant compared to that of g and that the effective 
mass may diverge at th e tra nsition and should increase 
with the magnetic field 258]. The strong increase of the 
effective mass near electron crystalliz ation also follows 
from Gutzwiller's variationa l me thod |259j |. which was 
applied to silicon MOSFETs |260|, and from the dynam- 
ical mean- field theory |26l| . Although the sharp increase 
of the mass is in agreement with the experimental re- 
sults, the suggested dependence of m on the degree of 
spin polarization is not confirmed by the data. 

Concluding this section, I would like to make some 
more remarks on the Fermi-liquid-related concepts. An 
idea was expressed to connect the observed effective mass 
enhancement to possible formation of a coupled-valley 
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state in bivalley electron systems |262t l263j | . Still, it is at 
odds with the fact that similar results are obtained for 
single-valley electron systems. An assumption was made 
that a plateau at the Fermi energy in the spectrum E(k) 
may form leading to a diverging effective mass (see, e.g., 
Ref. [264]). As for now, however, the dependence of the 
effective mass on temperature, resulting from the plateau 
formation, is not consistent with the experimental find- 
ings. A prediction that the electron density, at which the 
effecti ve m ass exhibits a sharp increase, is sensitive to dis- 
order 265] is not confirmed by the experimental data in 
available samples. After all, one can simply follow a clas- 
sical way of introducing phcnomcnologically Fermi-liquid 
parameters as the physical observables to be determined 
in experiment. 

IV. CONCLUSIONS 

A critical analysis of the available experimental data 
for 2D electron systems shows that consequences of the 
scaling theory of localization are not confirmed. The 
main points to be addressed by theory are the problem 
of finite bandwidth of the extended states in the Landau 
levels, and that of a true metal-insulator transition in 
zero magnetic field whose existence is strongly supported 
in low-disordered 2D electron systems, but remains un- 
certain in 2D electron systems with high disorder. Also, 
there is still no theoretical description of the oscillations 



of the metal-insulator phase boundary as a function of 
perpendicular magnetic field. 

In the past four years, significant progress was made in 
understanding the metallic state in strongly-interacting, 
low-disordered 2D electron systems. The state is remark- 
able by the strong metallic temperature dependence of 
the resistivity caused by electron-electron interaction ef- 
fects. The spin susceptibility measured in low-disordered 
silicon MOSFETs using different experimental methods 
shows a sharp increase and possible divergence at a finite 
sample-independent electron density n x close to the crit- 
ical density n c for the metal-insulator transition. This 
indicates that the metal-insulator transition in clean 2D 
systems is driven by interactions. Unlike the Stoner in- 
stability, the increase in the spin susceptibility is caused 
by the enhanced effective mass rather than the g factor. 
The effective mass does not depend on the degree of spin 
polarization, so the mass enhancement is not due to spin 
exchange. A similar increase in the spin susceptibility 
is observed in other 2D carrier systems. It remains to 
be seen whether or not it indicates the occurrence of a 
spontaneous spin polarization at a finite carrier density. 
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A. Widom for valuable discussions. The author is sup- 
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